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1. 


INTRODUCTION 


This  is  the  fifth  quarterly  report  documenting  the  work  performed  during  a  two  year  Phase 
n  STTR  activity  entitled  "Particulate  contaminant  formation  and  transport  in 
microelectronic  manufacturing  processes."  The  overall  objective  is  to  produce  a  charging, 
transport,  and  growth  simulation  (CTGS)  tool  that  can  be  used  effectively  by  equipment 
manufacturers  and  users  to  reduce  particle  generation  m  fabrication  systems.  The  project 
includes  collaborations  with  Prof.  Steven  Girshick  (University  of  Minnesota)  on  ion 
induced  nucleation  and  Prof.  Mark  J.  Kushner  (University  of  IHinois)  on  particle  charging 
and  transport. 

1.1  Project  Objectives 

The  work  to  be  performed  in  this  Phase  II  STTR  includes: 

1.  To  obtain  and  adapt  models  for  particle  charging,  particle  growth,  and  electrostatic 
and  ionic  drag  forces  on  particles  from  Prof.  Mark  J.  Kushner's  group  at  the 
University  of  Illinois  Urbana-Champaign  (UIUC). 

2.  To  obtain  and  adapt  models  for  ion-induced  nucleation  from  Prof.  Steven 
Girshick's  group  at  the  University  of  Minnesota  (UMN). 

3.  To  implement  the  above  models  into  CFDRC's  reactor  model,  CFD-ACE,  to  obtain 
solutions  for  local  particle  formation  and  transport  in  conjimction  with  other 
macroscopic  processes  such  as  fluid  flow,  heat/ mass  transfer,  plasma  transport 
and  chemistry. 

4.  To  perform  detailed  parametric  studies  to  assess  the  effects  of  various  operating 
conditions  on  particle  formation  and  transport. 

5.  To  demonstrate  the  integrated  model  in  collaboration  with  an  equipment 
manufacturer  to  reduce /eliminate  particle  effects  in  a  processing  system. 

6.  To  transfer  the  technology  to  the  industry  in  collaboration  with  SEMATECH  and 
its  member  companies. 

7.  To  prepare  a  final  report  documenting  the  work  performed  during  the  Phase  II 
study. 
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1*2  Project  Status 


The  charging,  transport,  and  growth  simulation  (CTGS)  tool's  graphical  user  interface 
(GUI)  has  imdergone  further  development.  The  portion  of  the  GUI  handling  the  charging 
and  transport  portions  of  the  tool  is  nearing  completion.  The  basic  tools  have  now  been 
established  with  future  work  focusing  on  making  it  commercial  grade.  The  moment 
model  used  for  the  growth  aspect  of  the  tool  has  undergone  further  development. 

1.3  Overview 


Section 2  describes  activity  in  GUI  development  and  GUI  solver  coupling.  Sections 
describes  progress  made  modeling  nucleation  pathways  for  use  with  the  moment  model 
as  reported  by  the  sub-contractor.  Section  4  discusses  the  future  plans  for  the  CTGS  tool 
development. 

2.  GUI  DEVELOPMENT  OF  PARTICLE  CHARGING  AND  TRANSPORT 

A  graphical  user  interface  (GUI)  for  the  CTGS  tool  is  undergoing  further  testing.  Some 
modifications  that  were  a  result  of  the  testing  include: 

•  The  toolbox  portion  of  the  GUI  is  updated  assuring  that  all  frames  controlled  by 
variables  are  packed  properly  when  reading  in  a  new  model  file. 

•  Data  entered  into  an  entry  box  is  checked  for  integrity.  Checks  are  made  for  type 
(real,  integer,  character),  as  well  as  maximum  and  minimum  values.  If  a  number 
entered  is  out  of  range  it  will  be  replaced  by  the  specified  maximum  or  minimum 
number.  If  no  maximum  or  minimum  for  a  variable  is  desired  it  can  be  set  to  nil 
which  will  set  no  limit  on  the  data  value  entered. 

•  Output  to  entry  boxes  is  now  formatted. 

•  The  progress  of  calculation  is  now  graphed  using  a  tk  canvas  widget.  The  graph 
tracks  particle  number  and  average  velocity  versus  simulation  time.  The  graph 
has  self  adjusting  axis  values. 

•  Version  control  and  a  time  stamp  were  added  to  GUI. 
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•  The  current  working  directory  is  now  changed  to  the  location  of  the  model  file 
(.nam)  so  fortran  output  files  get  written  to  the  same  location  as  the  model  file  as 
opposed  to  wherever  the  gui  was  initially  called. 

•  A  routine  has  been  written  to  handle  inputting  the  more  generic  PFG  file  format 
where  the  i,  j,  and  k  indice  do  not  necessarily  line  up  with  the  x,  y,  and  z-axis, 
respectively. 

•  The  application  can  now  be  stopped  before  the  scheduled  end  and  have  all  neces¬ 
sary  result  files  generated.  This  is  achieved  by  using  a  signal  file  which  the  main 
application  queries  every  iteration  loop. 

•  The  dual  roles  of  the  dtmaxstep  variable  were  separated.  The  variable  was  used 
to  control  the  burst  period  as  well  as  the  period  between  writes  to  the  SPY  file. 

•  The  inner  loop  timing  mechanism  was  rewritten  to  run  on  a  global  clock  instead 
of  relying  on  integer  loop  values.  Using  global  clock  allowed  a  reorganization  of 
timed  events  and  particle  introduction  during  the  simulation. 

3.  MOMMENT  MODEL 

The  following  is  a  report  written  at  the  University  of  Minnesota  documenting  work  done 
on  the  nucleation  term  to  be  included  m  the  moment  model  described  here  and  in  earlier 
reports. 

3.1  Introduction 

This  report  documents  the  progress  made  during  the  fifth  quarter  of  the  Phase  II  STTR 
study  titled  "Particulate  Contaminant  Formation  and  Transport  in  Microelectronic 
Manufacturing  Processes".  The  objective  of  this  project  is  to  transfer/adapt  particle 
nucleation  models  developed  at  the  University  of  Minnesota  and  incorporate  them  into 
the  macroscopic  reactor  flow  and  plasma  models  developed  by  CFDRC.  The  combined 
code  will  provide  a  comprehensive  computational  tool  to  simulate  various  interactions 
(between  gas  species,  charge  carriers,  and  dust  particles)  during  microelectronic 
fabrication  processes  and  to  facilitate  the  analysis/ optimization  of  industrial  reactors. 
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3.1.1  Need  for  Numerical  Modeling 

Particle  contamination  currently  occurring  in  microelectronics  fabrication  processes  has 
been  understood  and  controlled  largely  through  empirical  methodologies  and  post¬ 
process  observation/ analysis.  The  sources  of  particle  contamination  can  be  attributed  to 
feed  gases,  gas  phase  reactions,  and  solid  surfaces.  The  remedy  methods  have  been 
developed  through  hardware  experimentation.  Contaminant  formation,  growth,  and 
transport  in  neutral  and  charged  gas  environment  involve  many  complicated  phenomena 
such  as  nucleation,  coagulation,  diffusion,  convection,  condensation,  evaporation, 
ionization,  and  interactions  between  particles,  molecules,  atoms,  ions,  and  electrons,  as 
well  as  these  particles'  interaction  with  external  force  fields.  It  is  difficult  to  understand, 
through  hardware  experimentation  or  purely  theoretical  analysis,  the  effect  of  mutual 
interactions  between  these  phenomena  over  a  range  of  operating  conditions.  The  sizes  of 
the  particles  generated  in  plasma  systems  are  of  the  order  of  a  micron  or  smaller.  It  is 
difficult  to  observe  and  measure  the  formation  and  transport  of  these  particles  even  with 
the  aid  of  highly  sophisticated  visualization  equipment.  The  underlying  physics  and 
chemistry  governing  the  particles'  behavior  in  plasmas  have  not  been  completely 
imderstood  yet.  As  a  result,  it  is  prohibitively  expensive  and  time  consuming  to 
experiment  with  process  design  and  equipment  development.  In  this  regard,  numerical 
simulations  can  be  a  very  effective  tool  to  model,  in  a  systematic  manner,  the  influence  of 
various  operating  parameters  on  particle  formation,  growth,  and  transport  plasma 
systems.  Such  a  model  can  be  used  to  imderstand  the  many  complex  and  coupled 
interactions  in  plasmas,  to  interpret  experimental  measurements,  and  to  predict  the 
behavior  of  the  particles. 

Although  for  many  years,  some  fxmdamental  modeling  work  has  been  done  to  address 
issues  of  particle  formation  in  areas  such  as  air  pollution,  clean-room  contamination, 
particle  generation,  powder  production,  and  industrial  atomization.  In  most  studies, 
phenomena  such  as  nucleation,  transport,  conservation,  kinetics  and  discharges  are  often 
treated  separately.  No  attempts  have  been  made  to  integrate  a  particle  model  with  general 
purpose  reactor  models.  There  is  currently  a  need  for  coupling  particle,  plasma,  flow,  and 
chemistry  models  with  reactor  simulations  in  order  to  facilitate  reactor  design.  This 
approach  will  benefit  a  large  number  of  industrial  users.  In  order  to  simulate  particle 
behaviors  in  plasma  processing  over  a  range  of  operating  conditions,  a  numerical  model 
should  have  the  following  features: 

a.  the  ability  to  model  flow  in  complex  reactor  geometries  with  arbitrary 
boundary  conditions; 
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b.  submodels  for  gas-phase  and  plasma  chemistry;  and 


c.  submodels  for  aerosol  dynamics  including  particle  nucleation  and  growth. 

The  application  of  advanced  numerical  models  will  result  in  improved  micro¬ 
contamination  free  processes  and  novel  equipment  designs  for  semiconductor  device 
fabrication. 

3.2  Project  Objectives 

The  overall  objective  of  the  proposed  work  is  to  develop  and  transfer  advanced 
computational  codes  for  particle  nucleation  and  growth  from  the  University  of  Minnesota 
to  CFDRC.  These  models  will  be  integrated  into  an  existing  Computational  Fluid 
Dynamics  (CFD)  code  capable  of  solving  gas  phase  transport  and  chemistry  in  complex 
reactor  configurations.  A  low-pressure  plasma  reactor  model  has  been  preliminarily 
developed  at  CFDRC  and  is  subjected  to  refinement  and  completion  pending  further 
fimding.  The  combination  of  these  three  categories  of  codes  will  result  in  the  development 
of  a  overall  (comprehensive)  computational  tool  that  simulates  interactions  between  fluid 
flow,  discharges,  and  particulates  over  a  range  of  operating  conditions. 

Our  specific  goal  for  the  Phase  11  STTRis  to  develop  and  transfer  realistic  models  of  particle 
nucleation  and  growth  under  conditions  representative  of  microelectronic  manufacturing 
processes.  These  may  include  such  processes  as  homogeneous  nucleation,  ion-induced 
nucleation  and  chemical  nucleation,  i.e.  nucleation  resulting  from  chemical  reactions  of 
gas-phase  precursor  species. 

3.3  Project  Status 

The  following  tasks  have  been  accomplished  in  meeting  the  goals  of  the  Phase  II  STTR 
study:  An  atomistic  modeling  framework  has  been  developed  for  calculation  of  particle 
nucleation  rates  in  microelectronics  manufacturing  processes.  Classical  nucleation 

theory^*^  is  based  on  the  "capillarity"  approximation,  i.e.,  the  assumption  that  small 
clusters  are  "liquid"  droplets  whose  properties  equal  their  values  for  the  bulk  liquid,  and 
also  presupposes  the  existence  of  a  conderrsable  vapor.  However,  this  approach  is  not 
suitable  for  treating  nucleation  in  microelectronic  fabrication  processes,  where  the  critical 
cluster  size  is  often  very  small,  and  in  some  cases  the  nucleation  proceeds  via  chemical 
reactions  between  gas  phase  radical  species  in  the  absence  of  a  supersaturated  vapor.  In 
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the  atomistic  approach^’^,  nucleation  rates  are  calculated  based  on  the  properties  of  small 
discrete  clusters,  which  eventually  grow  to  large  aerosol  particles.  The  nucleation  rate 
expression  obtained  from  this  approach  is  in  the  form  of  a  summation  over  discrete  cluster 
properties.  The  summation  converges  rapidly  aroimd  the  critical  cluster  size,  which  is 
often  as  small  as  a  few  atoms.  Summation  expressions  may  be  derived  for  all  possible 
mechanisms  of  particle  nucleation  including  homogeneous  nucleation,  ion-induced 
nucleation,  and  chemical  nucleation.  This  approach  allows  the  direct  utilization  of 
experimental  and/or  computational  data  for  cluster  properties,  which  are  known  to  be 
poorly  approximated  by  the  traditional  "capillarity"  model.  Sample  calculations  have 
been  performed  for  the  nucleation  of  silicon  particles  in  microelectronics  fabrication 
processes  assuming  various  nucleation  scenarios  including  (1)  homogeneous  nucleation 
from  silicon  vapor;  (2)  condensation  of  neutral  silicon  vapor  onto  silicon  anions,  for 
conditions  representative  of  silicon  PECVD;  and  (3)  polymerization  of  silicon  hydride 
species  imder  conditions  representative  of  silicon  LPCVD  and  PECVD  processes.  The 
nucleation  rates  derived  from  these  atomistic  approaches  may  be  used  as  particle  source 
terms  in  moment  models  of  aerosol  dynamics  being  developed  by  CFDRC. 

3.4  Summary  of  Recent  Results 

In  past  quarters,  we  have  explored  several  clustering  pathways  (including  homogeneous 

nucleation  of  silicon  vapor,  ion-induced  nucleation  of  Si  anions,  and  chemical  nucleation 
routes)  for  the  formation  of  particles  in  LPCVD  and  PECVD  systems  used  for  depositing 
silicon  from  silane  precursor.  In  the  recent  past,  however,  there  has  been  considerable 
experimental  evidence  that  the  formation  of  particles  in  silicon  deposition  systems 
proceeds  through  gas-phase  chemical  reactions  involving  hydrogenated  silicon  cluster, 

i.e.  silicon  hydrides^'^^.  Several  reaction  sequences  involving  silicon  hydride  species  have 
been  proposed,  including  those  involving,  neutral  radicals,  and  even  cations.  Of  these  the 
reaction  pathway  appearing  to  have  the  strongest  experimental  support  is  the  anion 

clustering  route^®'^^  of  the  form: 


■5'„  +  SiH,  «  ^  (1) 

for  n=l,2,3... Nucleation  rates  for  this  reaction  sequence  have  been  previously  calculated 
using  the  atomistic  theory  for  hydrogenated  clusters  containing  up  to  5  silicon  atoms. 

Thermodynamic  property  data  for  these  clusters  are  now  available  in  the  literature^^. 


6 


4396/5 


In  this  report  we  present  results  from  alternative  chemical  nucleation  scenarios  involving 
neutral  species.  Two  possible  reaction  sequences  for  silicon  have  been  proposed  in  the 

literature,  and  could  apply  to  both  PECVD  as  well  as  LPCVD  processes®'^'^^: 

-  1^2n  +  ^^^2  ^  ^^n^2n  +  2  (2) 


-  1^2n  -  2  ^  ^^n^2n  ^2  (3) 

Following  the  usual  formulation  of  the  atomistic  nucleation  theory^  for  chemical 
nucleation,  nucleation  rates  may  be  calculated  using  expressions  of  the  form  (for  Reaction 
3), 


J  = 


^SiH2^SiHj\ 


M 

Z 

i  =  I 


y-1 


N 


exp 


H.J 


'Vj 


(4) 


Figure  1  shows  calculated  nucleation  rates  versus  temperature  for  concentrations 
reasonably  representative  of  silane  PFCVD  and  LPCVD  (NsiH2=  10^^  cm'^,  NgiH^  =  10^^ 
cm'3)  using  the  various  chemical  nucleation  model  expressions.  For  the  case  of  reactions 


(1)  and  (3),  the  "driving  force"  a  = 


SiH, 


N 


was  assmned  to  be  unity.  Approximate  rate 


constants  for  the  clustering  reactions  (1)  to  (3)  have  been  reported  in  the  literature^ ^42^ 
and  have  been  used  in  these  calculations. 


From  Figure  1  it  is  seen  that  at  lower  temperatures  (T  <  1000  K),  the  atomistic  rates  of 
chemical  nucleation  depend  strongly  on  the  reaction  sequence.  The  neutral  reaction 
sequences  yield  higher  nucleation  rates,  and  thus  appear  to  be  viable  even  in  PFCVD 
processes.  However,  it  must  be  noted  that  the  nucleation  rates  also  depend  on  the 
concentrations  of  the  "growth"  species  and  the  "initiating"  species  participating  m  these 
reactions.  At  present,  these  values  have  been  estimated,  and  may  be  inaccurate  by  several 
orders  of  magnitude.  In  order  to  draw  conclusions  regarding  the  most  likely  nucleation 
route,  it  is  necessary  to  couple  the  atomistic  nucleation  models  presented  here  to  gas-phase 
chemistry  models  being  developed  by  CFDRC. 
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Figure  3-1.  Calculated  Nucleation  Rates  of  Silicon  from  Chemical  Reaction  Pathways 

4.  FUTURE  WORK 


The  atomistic  approach  to  calculating  nucleation  rates  has  been  extended  to  clustering 
reaction  sequences  involving  neutral  hydrogenated  clusters  of  silicon.  The  next  step  is  to 
integrate  these  atomistic  expressions  into  aerosol  dynamics  models,  e.g.  a  moment 

model^^  capable  of  accounting  for  other  processes  such  as  particle  growth,  particle 
charging  and  neutralization,  particle  transport  and  deposition  onto  reactor  walls. 

The  current  status  of  the  project  is  shown  in  Table  4-1.  Future  work  will  address  both 

areas  of  the  technology  transfer  focusing  on  improvements  to  the  transport,  charging  and 
growth  models. 


e  moment  model  will  undergo  further  development  to  incorporate  the  atomistic 
nucleation  model  as  a  source  term.  The  possibility  of  forming  a  discrete-modal  model  will 
be  explored.  A  discrete-modal  model  solves  separate  General  Dynamic  Equations  (GDEs) 
for  small  clusters  and  the  moments  of  the  GDE  assuming  a  lognormal  size  distribution  for 

particles  above  a  specified  size.  Coupling  a  discrete  model  with  a  moment  model  will  be 
explored  as  an  option. 
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Table  4-1.  Status  of  Project  Tasks 


The  GUI  for  the  charging  and  transport  will  undergo  further  testing.  More  efficient 
information  sharing  from  ACE  and  PLASMA(ICP)  code  to  the  transport  code  will  be 
explored  such  as  through  the  restart  files.  The  transport  model  which  currently  assumes 
monodisperse  systems  will  be  altered  to  transport  more  than  one  particle  size  at  the  same 
time.  Other  modifications  which  have  yet  to  be  finalized  could  include  transport  of  non- 
spherical  particles  (i.e.,  flakes),  accoxmting  for  particle  coulomb-coulomb  interactions, 
making  particle  charging  time-dependent,  changing  the  charging  model  for  smaller  size 
particles,  and  allow  for  variance  in  size  as  particles  collide  with  each  other  increasing  in 
size  or  atoms  evaporate  off  of  the  particle  reducing  its  size. 
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